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We introduce a generalization of Obukhov's model [A.M. Obukhov, Adv. Geophys. 6, 113 (1959)] 
for the description of the joint position-velocity statistics of a single fluid particle in fully developed 
turbulence. In the presented model the velocity is assumed to undergo a continuous time random 
OO ' walk. This takes into account long time correlations. As a consequence the evolution equation 

for the joint position- velocity probability distribution is a Fokker-Planck equation with a fractional 
time derivative. We determine the solution of this equation in the form of an integral transform and 
derive a relation for arbitrary single time moments. Analytical solutions for the joint probability 
distribution and its moments are given. 
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I. INTRODUCTION 



The statistics of turbulent flows is described by phenomenological theories dating back to Kolmogorov and Obukhov 
[EliLIII]- In a Lagrangian treatment the path X(t, y) and the velocity U(£, y) of advected particles initially starting at 
the position y are determined by the acceleration A(£, y) which in principle is given by the Navier-Stokes equation: 

■ dU(t,y) . . 

\; y> — A(t, y) . 1 
>v dt 



The central statistical quantity is the joint position-velocity probability distribution of the particle: 

/(u,x,t)=<J(x-(X(i,y)-U(0,y)i-y))5(u-(U(i )y )-U(0,y))> , (2) 

^ ■ where the brackets denote a suitable average over a stationary statistical ensemble. Here and in the following we 
shall be interested in the statistics of increments u(t) = U(t, y) - U(0, y), x(t) = X(i, y) - U(0, t)t - y. The pdf © 
obeys the initial condition, /(u, x, 0) = 5(u)<5(x). Furthermore, Kolmogorov's theory of 1941 (K41) [l[ suggests the 
' existence of scaling behaviour 

/(u,x,t) = i/ (-^,A) • (3) 

Although such a scaling behaviour may not exist due to intermittency the assumption of normal scaling behaviour 
•i-h ' ([3]) usually serves as a first starting point. 

Due to the fact that a successfull statistical approach based on a treatment of the Navier-Stokes equation is still 
missing, one has to resort to heuristic assumptions about the statistics of A(t, y). An early model for /(u, x, t) has 

O been introduced by Obukhov 0J], who assumed, that the turbulent acceleration A(t,y) is a Gaussian, ^-correlated 

L" ' random force: 



X- <A l (t,y)A 3 (t',y)> = 2S l] S(t-t') . (4) 

H ' 

' In turn the probability distribution (J2J) obeys a Fokker-Planck equation: 

-/(u,x,i)+u-V x /(u,x,t)=A u /(u,x,f) . (5) 

The solution of the Obukhov-model with the initial condition /(u, x, 0) = <5(x)<5(u) possess a Gaussian form 
exhibiting scaling behaviour 

<u 2 {t)>~t , <x. 2 {t)>~t 3 , (6) 

which is consistent with the phenomenological scaling theory of Kolmogorov (K41). However Obukhov's model 
contradicts recent experimental results Q: the probability distribution for the velocity increment is far from 
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Gaussian and the increment behaviour can not be a simple random walk as suggested by Obukhov. The reason is 
the intermittent character of turbulent flows. Nonnormal statistics for Lagrangian variables may originate from long 
time correlations in the turbulent field. During the last few years lots of efforts have been put into the formulation 
of more sophisticated phenomenological theories which can take into account these facts (c.f. the review article 0]). 
Recently a connection between the velocity increment statistics of a Lagrangian particle and a type of continuous 
time random walk has been introduced [8j by a truncation of an infinite chain of evolution equations for multiple 
particle probability distributions. The following evolution equation for the single time joint position velocity increment 
probability distribution has been obtained 

Fl 1 rl t l dt' 

{_ + u-V x }/(u,x,f) = — - J _ V M Q(u 2 ) [V u /(u, x - u(t - f ), fl]u>=u , 

Q i3 -(u a ) = u^-^fe-Qr + Q^] . (7) 

(r(a) denotes the Gamma- function). In the case of isotropic turbulent flows, the matrix Qy (u ) has to be invariant 
with respect to rotations leading to the form given above. Furthermore, the variables Q\ and Qi are constant due 
to the fact that solutions of the above equation are required to allow for scaling behaviour of the velocity increment 
U ~ t 1 / 2 (for constant value of a). The parameter a is taken from the interval < a < 1. This equation generalizes 
Obukhov's model in several respects. First, it introduces a temporal memory. Second, the simple diffusion process 
in velocity space is changed to a diffusion process with velocity dependent diffusion coefficient. Third, retardation 
effects with respect to the spatial coordinate appear. It has been shown that solutions of the resulting equation for 
the velocity pdf G(u, t) = J dx/(u,x, t) yields excellent approximations to the experimentally determined velocity 
pdf's, provided the parameter a is allowed to vary with the time increment. A complete solution to the equation ([7]) 
has not yet been obtained. 

In the present paper we investigate a simpler phenomenological model by disregarding the retardation in the spatial 
coordinate: 



{^ + u-V x }/(u,x,i) = OJ D t 1 - Q i(u,V„)/(u,x,t) + ( 5(t) ( 5(x) ( 5(u) (S) 



Here, L(u,V u ) is a diffusion operator with a velocity dependent diffusion constant (compare eq. <[7j)). We have 
explicitly included the initial condition. The operator oD}~ a denotes the Riemann-Liouville fractional differential 
operator (see e.g. (l5j): 

This is a straightforward generalization of integer order differentiation to fractional orders. The presence of the 
fractional differential operator eq.® introduces temporal memory effects depending on the parameter a. We assume 
that < a < 1. In the limit a — * 1 our model eq.® reduces to the ordinary Obukhov model eq.© due to the 
property oD^F(t) = F(t) if one takes L(u, V„) to be the diffusion operator L = A u . 

The purpose of the present paper is not to state the accurate evolution equation for the Lagrangian position- velocity 
increment. This task has to be performed on the basis of a theoretical analysis of the Navier-Stokes equation, as has 
been started in [8[, or by sophisticated data analysis of the experimentally obtained Lagrangian path's of particles. 
An assessment of the underlying stochastic process, in our opinion, has to include the description of multiple time 
distributions. Such an analysis for a class of continuous time random walks has recently been started [9j. The purpose 
of the present paper is to consider an extension of Obukhov's model to a class of fractional diffusion equations, whose 
study quite seemingly is interesting by its own. It is hoped that the obtained results contribute to a more detailed 
understanding of Lagrangian turbulence statistics. 

The paper is outlined as follows. In the next section we focus on the probability distribution for the velocity 
increment and review how it can be determined by methods developed in the theory of continuous time random walks 
(for a recent review of the topic we refer the reader to [lj|, dU). In section III we shall show that the joint position 
velocity pdf can be expressed as an integral transform 

/(u,x,t) = J dz J dsh(z,x.;s,t)fi(u,z,s) , (10) 

where /i(u, x, s) is the solution of cq. ||5J) for a = 1 and h(z,x;s,t) a positive function. Special emphasis is put 
onto the case of L(u, V«) = A„. Here, the statistics of the velocity is characterized by subdiffusive behaviour related 
to temporal memory effects. It is shown that both velocity and position of the particle reveal anomalous diffusive 
properties. We explicitly determine the probability distributions of velocity and spatial coordinate. For a diffusion 
operator L(u, V u ) which allows for scaling solutions /(u, x, s) we discuss how normal scaling behaviour < u 2 (t) >^ t 
may arise in connection with nonnormal statistics. 
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II. THE PROBABILITY DISTRIBUTION OF THE VELOCITY INCREMENT 

After integration over the position variable in eq. ([8]) an evolution equation for the marginal probability distribution 
of the velocity G(u, t) = J tfx/(u, x, t) is obtained in the form of the fractional diffusion equation: 

J^G(u, t) = Di~ a L(u, V u )G(u, t) + S(t)S(u) . (11) 

If we take L(u, V„) = A u we obtain the fractional diffusion equation for a continuous time random walk [l(|, [ill ]. 
Accordingly the second order moment reveals subdiffusive characteristics: < u 2 (t) >~ t a . The corresponding pdf is 
well-known and will be given below. 

The solution of the fractional diffusion equation pip is conveniently expressed as an integral transform [l2| , [l3| 

poo 

G(u,t) = dsh a (s,t)Gi(u,s) , (12) 
Jo 

in which the function G\ (u, s) is obtained as solution of the ordinary diffusion equation (|11[) with a = 1 and the 
integral kernel h a (s,t) is given as the single time probability distribution of an inverse Levy-stable process [lH, [l3j : 

ft »(*,*) = i™^af4 J ) ■ (13) 



Here L a (t) denotes a one-sided Levy-stable distribution of order a. This may be shown by inserting the ansatz (|12p 
into eq. (fTT|) leading to the solvability conditions 

g- t h a {s,t) = - a Dl~ a — h a (s,t) 
oDl- a h a (0,t) = S(t) . (14) 

These two equations determine h a (s,t) in the form eq. ([13]) . The properties of the diffusion operator L(u, V u ) enter 
via the pdf Gi(u, t). The form of the pdf, eq. (fl"2"]) has the following interpretation. There is a usual diffusion process 
w(s) described by a Fokker-Planck operator L(w, V m ) with respect to an intrinsic time s. Additionally, there is a 
prozess s = s(t) relating physical time t and intrinsic time s with a probability distribution h a (s,t). The random 
process u(i) is given by u(t) — w(s(f)). 

III. SOLUTION OF THE GENERALIZED OBUKHOV-MODEL 

Motivated by the results for the fractional diffusion equation we look for solutions of eq. ((SJ) which can be expressed 
as an integral transform similar to eq. (|12|) : 



/oo />oc 
dz / ds/i(z,x;s,t)/i(u,z, s) . (15) 
-oo JO 

/i(u, z, s) denotes the solution of the ordinary Obukhov- model eq.® specified by a = 1: 

d 

— /i(u,x,s) +u- V x /i(u,x,s) = L(u, V„)/i(u,x,s) . (16) 

The structure of the pdf, eq. ([T5]) . has the following interpretation. Let us consider the case L = A u . There is a 
stochastic process with respect to the intrinsic time s given by 

dz 

T = u ( s ) 
ds 

dw , , . 

— - F(w,,s) . (17) 

Furthermore, there is a process s(t) relating physical time t and intrinsic time s. The velocity u(i), thereby, is given 
by u(i) = w(s(t)). If the space variable x(t) would be obtained via x(t) = z(s(t)) the function h would not depend 
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on x and z and would be simply given by the function h a (s, t) of the preceding section. However, the space variable 
x(t) is defined by the relationship 



<2x(£) 
dt 



u(t) 



(18) 



This explains the fact why /i(z, x; s, t) depends on the two spatial variables x, z but not on the velocity u. 

In the following we shall derive a fractional differential equation for h(z,x; s, t), corresponding boundary conditions 
and obtain an explicit solution. We proceed as follows: the ansatz ea.([T5|) is substituted into the evolution equation 
([8]). The operator L(u, V„) only acts on /i(u, z, s) such that eq. (Tl6|) can be applied: 



dz 



ds 



d_ 

at 



u ■ V x h(z, x; s, t)fi (u, z, s) 



J dzj^ dsh(z,x;s,t) f^+ u -Vz) /i(u,z,s)+<J(i)<J(u)«J(x). 



(19) 



In the integral on the right hand side we perform a partial integration with respect to z and s. Since /i(u, z, s) 
vanishes for |z| — > oo the following boundary terms are obtained: 



dzh(z, x, s, t)/i(u, z, s) 



(20) 



s=0 



We can assume /i(z, x; s = oo, <) = 0. 

Let us consider the boundary term using the initial condition /i(u, z, s = 0) = 6(u)S(z). In that case, the boundary 
term ^ yields: 



dzh(z, x; s, i)/i(u, z, s) 



-oA 1 " Q ^(0,x;0,f)<5(u) = -5(x)J(u)<y(t) 



(21) 



s=0 



This term cancels the last term on the right hand side of (fT9|) provided we postulate the validity of the boundary 
condition 



D\- a h(z = 0, x; s = 0, t) = 5(t)6{x) 
As a consequence, the function /i(z,x;s, t) has to obey the following equation 



(22) 



d 



+ u • V x I h(z, x; s, t) = - Q D] a ( ^-+u- V z ) h{z,x;s,t) 



dt 

a Dl- a h(z = 0,x;s = 0,t) = S(t)S(x) 







ds 



(23) 



This is the generalization of eq. (TT4"|) . 

We have to add the following remarks. First, equation (|2"3"|) determines the u-independent function h(z, x;s,t) 
by a relation including the velocity u. The fact that there is no u-dependence has been crucial for arriving at the 
equations (|23|) . Below, we shall show that we actually can find a u-independent solution. Second, we point out that 
the properties of the diffusion operator L(u, V u ) do not show up in the determination of the function h(z, x; s, t). The 
properties of this operator are included in the pdf /(u, z, s) and we have obtained the solution in terms of an integral 
transform for a large class of stochastic processes. 



A. Determination of h(z, x; s, t) 

The solution of differential equations containing the Riemann-Liouville fractional differential operator is simplified 
by changing to Laplace-space. This is due to the fact that the integral in the definition eq.© is actually a Laplace- 
convolution. Throughout this discussion we denote Laplace-transforms as follows: F(X) := £{F(t)} := J °° dte~ xt F{t). 
By performing a Laplace-transformation of eq. (|23p we derive the following first order partial differential equation: 

(A + u • V x ) h(z, x; s, A) — h(z, x; s, 0) = - A 1 "" ( ^- + u • V z ) h(z, x; s, A) 



.s 



+ (^ + u.V z )[oA 1 ^M Zl x; S ,t)] t=0 



(24) 
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The last term on the right hand side is due to the partial time derivative in the fractional differential operator. It can 
be set to zero: [oD^ a h(z, x; s, t)] . = 0. Equation (pM|) is valid for s > 0. For s = 0we have the boundary condition 

h(z, x; s, 0) = S(z — x)S(s). Thus we obtain the final form of the evolution equation for h(z, x; s, A): 

^- + u-V a + A Q - 1 u-V x ih(z,x;s,A) = -X a h{z,x; s, X) . (25) 



Linear partial differential equations of this type can be solved by the method of characteristics (see e.g. [16j). This 
method will be applied below. 

If we choose s as parameter of the characteristics, eq. ([25l) can be written as 



^(z( S ),x( S );A, S ) = |^ + ^.V z + ^-V x U(z( S ) )X ( S );A lS ) 

= -A" h(z(s),x(s);X,s) . (26) 

Clearly the solution of this differential equation reads 

h(z(s),x(s);X,s) = h(z(0),x(0); A,0) e'^ 3 . (27) 

The characteristics z(s), x(s) are determined by ordinary first order differential equations: 

dz(s) . 

— : — = u — ► z(s) = us + z , 
as 

A"" 1 !! -> x(s) = A a_1 us+x . (28) 



(is 



For the initial condition h(z(0), x(0); 0, A) = /i(zo, xo; 0, A) we specify a sharp distribution at xq. The dependence 
on A is chosen as A Q_1 to be consistent with the boundary condition ([23)1 . If we put Zo = the inital condition takes 
the form: 

h(zo = 0,x o ;0,A) = (5(x )A Q - 1 = S(x- X a - 1 z)X a ~ 1 . (29) 



Finally we obtain the kernel of the transformation eq. (|T5|) as solution of the partial differential equation (f25|) in 
Laplace-space: 

/i(z,x;s,A) = «J(x- A Q - 1 z)A Q - 1 e- A ° s = 5(x - A Q " 1 z)/i Q (s, A) . (30) 

In the last step we substitute the Laplace-transform of h a (s,t): £{h a (s,t)} — A a_1 e _A ° s . 

From this expression the integral kernel in physical time h(z, x; s, t) is given as inverse Laplace-transform: 

h(z,x;s,t)=£- 1 {5(x-X a - 1 z)}*h a (s,t) . (31) 

The asterisk * denotes a Laplace-convolution with respect to t. However due to the A-dependence of the (^-function a 
closed form of the inverse transformation £ _1 {i5(x — A Q_1 z)} could not be calculated and the explicit expression of 
h(z, x; s, t) remains unknown. 

Still the transformation eq. ([T5|) provides a useful equation for characterizing the solution /(u, x, t). Since inverse 
Laplace-transform and integration commute we can give an expression for the transformed solution f(u, x, A) := 
£{/(u )X) i)}: 

/oo />oo 
dz / ds h(z, x; s, A)/i(u, z, s) . (32) 
-oo Jo 

Thus statistical quantities can be calculated in Laplace-space and then transformed to physical time. It will become 
clear in the following that with the help of the integral kernel h(z,x; s, A) the solution of the generalized Obukhov- 
model can be completely characterized. 

However, several remarks are in order. First of all, one can easily check that the solutions /(u, x, t) of eq.© are 
normalized: 

/oo />oo />oo 

dx I duf(u,x,t)= du/(u,t) = l . (33) 
- oo J — oo J — oo 
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Additionally one states that h(z, x; s, A) is always nonegative and so are h(z, x; s, t), /i(u, z, s), and, in turn, /(u, x, t). 
This proofs that the generalized Oboukhov model actually defines a probability distribution for the case < a < 1. 

Second, we have to show that the above assumptions are fullfilled. To this end we note that the quantity h(z, x; s, t) 
in fact is independent on the velocity u. Furthermore, we have the relationship 

h(0, x; s, A) = 5(x)h a (s, A) , (34) 

which in real space reads 

Dt- a h(0,x;0,t) = 5(x) Dt- a h a (Q,t) = S(x) . (35) 

Since h a (s, 0) = <5(s) we obtain 

h(z,x;s,0)=£- 1 {5(x-X a - 1 z)}*S(s)~6(s) . (36) 

The case a = 1 has to lead to the ordinary Obukhov model. Evidently the integral kernel eq.(|30p in this limit is given 
by the Laplace transform 

h(z,x;s,X) = <5(x - z)e- As . (37) 
The inverse Laplace-transformation is easily calculated: 

h{z,x;s,t) = S(x-z)S(t- s) , (38) 
such that the transformation cq.([T5]) leads to the solution of the ordinary Obukhov- model. 

The second limit case is supposed to yield the solution G(u, t) of the fractional diffusion equation. Perform- 
ing the integration of eq. ([3"2")) with respect to x results in: 

/OC />OC pOO pOO 

dx dz ds6{x- A"- 1 z)/i Q (s, A)/i(u,z,s) = / ds h a (s, X)Gi(u, s) . (39) 
-oo J -oc JO JO 

The right hand side of this equation is the Laplace-transform of the solution eq. IT2|) of the fractional diffusion equation 
(fTTj) . The agreement in this case is achieved due to the choice A Q_1 in the initial condition eq. (p?9"l) . 



B. The joint probability distribution /(u, x, t) 

In the following we derive a formal expression for the joint probability distribution h(z,x; s, A). First the Laplace- 
transform of /(u, x, t) is calculated using the transformation eq. 



/oo />oo 
dz / ds 5(x — \ a ~ 1 z)h a (s, A)/i(u, z, s) 
-oo JO 

ds A d(1 ~ Q) /i Q (s,A)/i(u,A 1 - Q x,s) . (40) 



Here the integration with respect to z is performed in d dimensions. Let us assume that the diffusion operator 
L(u, V M ) allows for solutions exhibiting scaling behaviour of the form 

/i(u,z,a) = s(1 l 2S)d H(^,^) (41) 

This yields: 

r°° 1 \l-a u 

/(u,x,A)=/ dsX d ^h a{sA) ^_^^ H{ } . (42) 



/ • • • s (l+2<5)a • s it» s u 

Let us consider the case of the simple diffusion operator L = A. /i(u, z, s) is then explicitly given as: 

s ( V3 Y f u 2 3z-u 3z 2 l 
/1(U ' Z ' S) =( V 2^J ^{-- + — ~-) • (43) 

This pdf exhibits scaling behaviour of the form (|4"Tj) with 5 = 1/2. 

We have not determined the Laplace inversion of this formula. However, we shall obtain the probability distributions 
for the velocity u and the spatial variable x. 
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C. The probability distribution G(x, t) 

We shall start with the probability distribution G(u, t) for the velocity u which is obtained by an integration with 
respect to the spatial variable x: 

G(u,t) = J dx/(u,x,i) . (44) 

In Laplace space, we obtain 

G(u,A)= / dsX^e-^'Gifaa) , (45) 



which, by Laplace inversion, yields 



/>oo 

G(u, t) = / dsh a (s, f)Gi(u, s) 
Jo 



(46) 



If we assume that the function G\ (u, s) exhibits scaling behaviour, 



= isG) > ( 47 ) 



r s 

the result may be reexpressed in a form which just evidences this scaling behaviour 



G(u,i) = / dah a (a)^g(^ s ) = ^g(^) . (48) 



Let us consider the case for L = A u . Here, S = I 



2 ' 



G(u,t) = Ids-— i—Lai^) J— . (49) 

This may also be represented according to 

r (j da / 2 u 2 a 

G{*,t) = Jd*L a { f ,)y=^e-*'> . (50) 

Formula ([48)) yields an interesting structure. First of all, it represents a probability distribution as a superposition of 
pdf's Gi(x, s) with different variances. Such a representation has been introduced by Castain et al. [l2j in a description 
of the intermittent behaviour of the Eulerian velocity increment pdf. Second, scaling behaviour < u 2 (t) >~ t is quite 
often related to a random walk behaviour in velocity space. However, the formula (|48p shows that such scaling 
behaviour may also arise for the cases Q where the scaling indices a, 5 are related by 

S<*=\ ■ (51) 
This e.g. happens for the rotationally symmetric solutions of the diffusion operator 

L(u,V u ) = V(u 2 )( 1 ^V . (52) 

Here, the pdf G(u, s) is given by 

G(u,s) = N{a)-^e~^ (53) 

These solutions have been shown to yield accurate reconstructions of the velocity increment pdf for turbulent flows 
measured in [B[ by adjusting the parameter a Q. For the case of turbulent flows, the determination of the parameter 
a and especially its time dependence is an open problem. 
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D. The probability distribution F(x, t) 

Let us now turn to the probability distribution F(x,t) for the position x. This distribution is obtained from the 
joint pdf /(u, x, t) by integration with respect to the velocity variable: 

F(x,i) = J du/(u,x,t) . (54) 

The Laplace transform of this quantity can be determined from (|40[) 

/>oo 

F(x,A)=/ ds\ d( > 1 - a ')h a {s,\)F 1 (\ 1 - a x,s) . (55) 
Jo 

If we consider a system exhibiting scaling behaviour, 

M*,s) = ^)F(-^) , (56) 

we may perform the transformation 

S = a \-V-«)/<l+S) (57) 

leading to 

F(x,A) = ^dSA^^^M^^/^A)^^^ (^_) . (58) 
Laplace inversion yields the probability distribution 

r°° l x 

The decisive step thereby, is the possibility to calculate the Laplace inverse of the function 

£ -x {A (i- a )/(i +5 ) /ia(s , A ( 1 - a )/( 1+5))A)} = £ - 1{A ^i-i e -^s } 

= V.^ : i t i+s L ^±i(— r+r) = hs«±x(s,t) . (60) 

ad + 1 S !+^TT !+« S S3TT 5 + x 

We can now explicitly state the result for the diffusion operator L = A u : 

One observes that the limiting case a = 1 leads to pdf of the Obukhov model. Furthermore, it is evident that the 
probability distribution of the spatial variable x is given by an expression similar to the one obtained for the velocity 
variable in form of an integral transform: 



Ffrt) = J daL [2+ay3 {a)^^^fa^/\^^ a . (62) 

This expression has to be compared with the one obtained for the pdf of the velocity ([50)1 . 

We mention that the solution (|58[) can be related to a continuous time random walk for the position increment x(t) 
of a Lagrangian particle. Continuous time random walks for the description of the distance between two Lagrangian 
variables in connection with Richardson's law have been discussed by Shlesinger et al. [l7l ]. 



E. Moments of the generalized Obukhov-model 

Starting from eq. (|40p we can derive a general relation for moments of arbitrary order of the joint probability 
distribution /(u, x, t). 

Single-time moments < a;"(t)it™(t) > are defined as follows: 

/oo />oo 
du dxi>™/(u,x,i) . (63) 
-OO J —OG 
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Therefore we obtain in Laplace-space with the help of eg. ([32)1 and eg. ([30)) : 

/oo />oo />oo 

du dz dsS(x- A Q - 1 z)/i a (s,A)x>™/i(u,z, s ) 
-oo J — oo J 



oc 







^(a-l)n / ds h a ( S: \) du dzzX/i( U ,Z,s) . (64) 



In the last step the integration over the 5-function has been performed. Here we recognize that simply the moments 
of the ordinary Obukhov-model appear: J du J dz u™z™ /i(u, z, s) — < z™(s)m™(s) >. As as consequence the inverse 
Laplace-transform leads to a general expression for the moments of the generalized model: 

/•OO 

<x?(t)u™(t) > = D { t a - 1)n dsh a (s,t) <z?(s)u?(s)> . (65) 



Jo 

Thus the moments < xf (t)u" 1 (t) > can be determined from the moments of the ordinary model < z™(s)u™(s) > by 
an inverse Levy-transform and fractional integration with respect to t. Since < zf(s)uV rl (s) > take the form of simple 
polynomials in s this results in fractional integration of single-time moments of the inverse Levy-stable process Q . 
These calculations can always be performed in an elementary manner. 

For example, by applying eq. (|65|) to the second order moments for the case of L = A u we obtain: 

< u 2 (,s) > = 2s - < u 2 (t) > = - - ' t a , 

1 (a + 1) 

< z 2 (*) > = |* 3 - < x 2 (t) > = T^rX+i , 

S 1 (a + 3) 

< z( S )u(s) > = s 2 - <x(t)u(t)>= * J a +\ (66) 

1 (a + 2) 

Clearly, in each case the limit a — > 1 of the Obukhov-model is satisfied. Therefore the generalized Obukhov-model 
eq. © with L = A„ is characterized by subdiffusive behaviour of the particle velocity corresponding to the fractional 
diffusion equation (fTI]) and mean square displacement < x 2 (i) > ~ t a+2 of the particle position. 



IV. CONCLUSION AND OUTLOOK 



We have introduced a non-Markovian model for the joint position-velocity probability distribution of a random 
walk particle as a generalization of Obukhov's model. Long-time correlations of the particle velocity are taken into 
account by fractional time derivatives. As in the limit case of the fractional diffusion equation, a solution can be 
found as an integral transformation of the Markovian solution. The integral kernel has been obtained by solving a 
fractional differential equation and yields a complete characterization of the single time statistical properties. As final 
results we have given analytical expressions for the joint probability distribution and its moments. The generalized 
Obukhov-model may thus serve as an appropriate model for random walk processes revealing anomalous diffusion in 
velocity- and position-space. 
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